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(57) Abstract: A method of optimizing a layout of selected parts to be cut from a piece of wood uses data representing geometric 
and defect-related characteristics piece and data representing geometric and grade characteristics. A subset of parts characterized 
by predetermined grade value and associated with a predetennined group of dimension values is defined, as weU as an arrangement 
of subdivided piece surface sections to be obtained through primary cutting operation. Subsections included in each section are 
defined according to geometric and defect-related characteristics, and a plurality of arrangements selected from the subset of parts 
to be obtained through secondary cutting operation, arc defined. Yield values associated with the arrangements arc estimated and 
compared to select the arrangement having a highest yield value. A basic yield value for the arrangement is estimated, and the 
optimization sequence repeated for new arrangements to estimate basic yield values, to select the best arrangement. 
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METHOD OF OPTIMIZING A LAYOUT OF SELECTED PARTS TO BE CUT 

Field of the invention 

The present invention relates generally to the freld of product 
5 manufacturing optimization, and more particularly to methods for optimizing a 
layout of selected parts to be cut from pieces of raw material, such as employed 
in the lumber processing industry. 
Brief description of the background art 

Product manufacturing optimization techniques have been developed 

10 during the past years to improve productivity associated with industrial processes 
such as manual and automated part cutting or sawing. Generally, such 
optimization techniques allow the detemiination of the arrangement of parts that 
provides the highest yield in terms of raw material usage and/or economic value 
from a given piece of raw nlaterial characterized by the presence of defects that 

15 have previously been qualified and located by spatial coordinates within a 
reference system, which can be either bi-dimensional or tri-dimensional to 
consider the two opposed main faces of the piece of raw material or the two pairs 
of opposed surfaces, depending upon the specific application considered. 
Typically, information about a set of parts to be cut from a piece of raw material is 

20 defined by a processing bill or order, which information includes requested 
quantities for each type of parts, quality grade associated to each requested 
product, geometrical characteristics such as length and width, unit cost, 
identification code, etc. Product cutting optimization is usually performed by a 
computer software which seeks to comply with a requested quantity of parts as 

25 set by the cutting bill, while optimizing the spatial arrangement of requested parts 
in order to maximize raw material usage and/or economic value obtained from 
the piece of raw material employed. Typically, an optimization software may 
provide production information details such as raw material yield, economic value 
yield or pricing, cutting layout, parts distribution, required number of cut 

30 operations, etc. Known optimization computer software is typically designed to 
be integrated into production equipment such as automated sawmill machinery or 
to be used as stand-alone systems for performing simulation from production 
data, in lumber processing, especially regarding hardwood cutting into wooden 
products such furniture components, panels and flooring pieces, crosscut-first 

35 rough mill operation has been considered as tiie main type of cutting process. 
However, shortage in large pieces of timber, raw material pricing increase and 
trends in the lumber industry toward diversified products as well as higher 
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production rates have oriented some sawing mills toward rip-first cutting 
processes, so that both crosscut-first and rip-first processes are now basic 
sawing processes employed in existing sawing mills. Crosscut-first and rip-first 
processes having their respective advantages and drawbacks, and depending 
5 upon the production context within which they are applied, product cutting 
performance observed for each of these main processing approaches may vary, 
rendering the task of cutting process selection a quite difficult one, as extensively 
discussed in the literature, such by Wiedenbeck, J.K. in "Dedding between 
Crosscut and Rip-first processing". Wood and Wood Products, August, 2001. 

10 Generally, rip-first process provides efficient elimination of defect areas adjacent 
the edges of the raw material piece, being more easily eliminated within a single 
section, while cross-cut process provides effident elimination of defect areas that 
are located over a short length portion of the piece by removing the whole 
defective portion. Most of the usual sawing processes favours the segmentation 

15 of each piece into a plurality of sections, such as rip-cut or crosscut sections, 
optimizing raw material surface according to a first reference axis, as discussed 
by Thomas, R.E, in "ROMI-RIP version 2.0: a new analysis tool for rip-first rough 
mill operations". Forest Products Journal, vol. 49, no.5, p. 35-40, 1999. Such 
sections are then further processed to eliminate defect areas while complying 

20 with dimensional requirements, such as discussed by Thomas, R.E. et al in 
"Dedsion-support software for optimizing rip-first and chop-first systems" Scan 
Pro, 8*** international conference on scanning technology and process 
optimization for the wood products industry. Flooring wood processing, which is 
a particular type of rip-first process wherein length dimension is undefined in 

25 cutting bills, and panel processing, which is a particular type of cross-first 
process wherein width dimension is undefined, bo\h offer advantages relevant to 
one-dimensional freedom axis for arranging parts, such one-axis optimization 
technique is disclosed in U.S. Patent no. 4.221,974 entitled "Lumber inspection 
and optimization system" issued September 9, 1990 to Muller et al. Another one- 

30 axis optimization method is disdosed in U.S. Patent 4,163,321 entitled "Method 
for sequencing the cutting of elongated stock" issued to Cunningham on August 
7, 1979, There is disdosed a method for cutting optimization of elongated stock 
such as boards of lumber having random unusable lengths containing defects, 
which method involves products order requirements, waste factors detemnination, 

35 probability assessment of having a given length of usable stock in each grade 
being processed, and from the infonnation above, determination of the 
arrangement of parts to be cut at one time which results in the lowest waste for 
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the entire cutting. Another known one-axis optimization method is disclosed in 
U.S. Patent 4,01 7,976 entitled "Apparatus and method for maximum utilization of 
elongated stock" issued to Ban- and al on April 19, 1977. which employs a yield 
optimization approach for crosscutting of usable lengths of stock such as boards 
5 of lumber having random unusable lengths containing defects, wherein lengths of 
stock required are determined, information describing the required lengths is 
stored, a value factor for each required length is calculated and stored ds well as 
statistical data describing usable lengths in various grades of stock, infomnation 
on various grades proportions being processed is stored, priority factors to 

10 increase the probability of cutting desired lengths are determined and stored, and 
finally, based on the infonrnation above, the positions of the backgauge indicators 
are calculated, which represent optimum possible combinations of lengths for 
each usable length which can be cut into the desired lengths, which position 
indicators are printed in full scale on a continuous sheet of paper. 

15 In order to improve the performance of optimization over known one-axis 

techniques, two-axis optimization methods and software have been developed, 
such as rip-first software discussed by Thomas, R.E. in "ROMI-RIP version 2.0: a 
new analysis tool for rip-first rough mill operations". Forest Products Journal, vol. 
49, no.5, p. 35-40, 1999, as well as crosscut-first software also discussed by 

20 Thomas in "ROMI-CROSS: An analysis tool for crosscut-first rough mill 
operations". Forest Products Journal, Vol. 48, no. 3, pp. 68-72. Such optimization 
systems allow the generation of an optimal cutting solution which consider two 

reference optimization axis simultaneously, which solution includes data on 

i 

selected anangements of parts, cutting position as well as output yield. ROMI- 
25 RIP and ROMKJROSS are two-axis optimization software products performing 
successive optimization steps considering quality/grades of products as defined 
in the cutting bill, from the highest grade to the lowest grade. A similar two^axis 
optimization technique is disclosed in U.S. Patent 3,329,181 entitled "Apparatus 
and method for cutting assorted lengths from material having irregular and 
30 random defects" issued July 4, 1967 to Buss and al. Another two-axis 
optimization technique is disclosed in U.S. Patent 4,805,679 issued to Czinner on 
February 21, 1989. which makes use of an expert system for performing the 
optimization task. Another two-axis apparatus for optimizing the yield of usable 
piece from boards and the like is disclosed in U.S. Patent 3,942,021 issued to 
35 Barr and al on March 2, 1976, which adopts a yield optimization approach 
wherein a primary bit matrix corresponding to a pattem of scanned unusable 
defect containing areas of a processed board of lumber, is formed and stored in 
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computer. Predetermined combinable unusable defect containing areas as well 
as predetennined combinable unusable non-defect containing areas are 
identified from the primary bit matrix, and the identified combinable unusable 
defect and non-defect containing areas are merged to produce a list defining a 
5 pattern characterizing usable areas. Predetermined billing requirements are 
established and stored, various cutting patterns for the workpiece are 
successively determined on the basis of usable area information and billing 
requirements, and the cutting pattern producing the optimum yield for a 
workpiece is selected. 

10 Although considered as a promising alternative to conventional one-^s 

optimization systems in use by rough-mill operations, most of existing two^s 
optimization systems are conceived on the basis of optimization methods that do 
not provide sufficient flexibility to ensure maximum optimization performance in 
specific applications, and therefore, they have not yet achieved general 

1 5 acceptance within the lumber industry. 
Brief summary of the tnvention 

It is therefore a main object of the present invention to provide a method 
of optimizing a layout of selected parts to be cut which allows flexibility of use 
while ensuring reliable optimization results. 

20 According to tiie above object, from a broad aspect, there is provided a 

method of optimizing a layout of selected parts to be cut from a piece of raw 
material witii respect to first and second orthogonal reference axis, using data 
representing geometric and defect-related characteristics of said piece and data 
representing geometric and grade characteristics of a set of parts to be cut The 

25 method comprises the steps of: i) defining a subset of the set of parts 
characterized by a predetermined grade value and associated with a 
predetennined group of first dimension values defined with respect to the first 
axis; ii) generating data defining an arrangement of subdivided piece surface 
sections to be obtained through a primary cut operation with respect to the 

30 second reference axis and according to one or more of the first dimension 
values; iii) generating data defining one or more subsections included in each 
piece surface section according to the geometric and defect-related 
characteristics of the piece; iv) generating data defining for each subsection a 
plurality of an^ngements of parts to be included tiierein and selected flrom tiie 

35 subset of parts, to be obtained ttirough a secondary cutting operation witii 
respect to the first reference axis; v) estimating yield values associated with the 
arrangements of parts; vi) comparing the yield values to select one of the 
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arrangements of parts having a highest yield value; vii) estimating a basic yield 
value the arrangement of subdivided piece surface sections; viii) repeating steps 
ii) to vii) for one or more new arrangements of subdivided piece surface sections 
to estimate corresponding basic yield values; and ix) comparing all basic yield 
5 values one witii anotiier to select the arrangement of subdivided piece surface 
sections associated with the arrangements of parts providing a maximal basic 
yield value to be included in the optimized layout of selected parts to be cut. 

According to the same main object, from a further broad aspect of the 
invention, there is provided a software product data recording medium in which 

10 program code is stored, which program code will cause a computer to perform a 
method of optimizing a layout of selected parts to be cut from a piece of raw 
material with respect to first and second orthogonal reference axis, using data 
representing geometric and defect-related characteristics of said piece and data 
representing geometric and grade characteristics of a set of parts to be cut, 

15 wherein tiie method comprises Uie steps of: i) defining a subset of the set of 
parts characterized by a predetermined grade value and associated witii a 
predetenmined group of first dimension values defined with respect to the first 
axis; ii) generating data defining an arrangement of subdivided piece surface 
sections to be obtained through a primary cut operation with respect to the 

20 second reference axis and according to one or more of tiie first dimension 
values; iii) generating data defining one or more subsections included in each 
piece surface section according to the geometric and defect-related 
characteristics of the piece; iv) generating data defining for each subsection a 
plurality of arrangements of parts to be included therein and selected from ttie 

25 subset of parts, to be obtained through a secondary cutting operation with 
respect to the first reference a)ds; v) estimating yield values associated with ttie 
arrangements of parts; vi) comparing the yield values to select one of tiie 
arrangements of parts having a highest yield value; vii) estimating a basic yield 
value the arrangement of subdivided piece surface sections; viii) repeating steps 

30 ii) to vii) for one or more new arrangements of subdivided piece surface sections 
to estimate conresponding basic yield values; and ix) comparing all basic yield 
values one witii another to select the arrangement of subdivided piece surface 
sections associated with the arrangements of parts providing a maximal basic 
yield value to be included in tiie optimized layout of selected parts to be cut. 



wo 2004/051523 



PCT/CA2003/001872 

6 



Brief description of the drawings 

A preferred embodiment of a method of optimizing a layout of selected 
parts according to the present invention will now be described in detail with 
respect to the accompanying drawings in which: 
5 Fig. la is a process flow diagram representing the main method steps as 

part of the prefenred embodiment of the method according to the invention; 

Fig. 1b is a process flow diagram complemen^ry to the diagram of Fig. 
1a, representing further method steps involving a group of n predetermined 
grade values of decreasing cutting priority order and a further optimization mode 
10 to be considered to select the arrangement of parts providing a maximum basic 
yield value; 

Fig. 1c is a process flow diagram complementary to the diagram of Fig. 

lb, representing a sequence of steps associated v^th a recycling function of the 

method according to tiie preferred embodiment, according to which one or more 
15 additional parts to be included in free areas of a piece sub-section are selected 

from a furtiier subset of parts; 

Fig. Id is a process flow diagram complementary to the diagram of Fig. 

1c, representing a case where a further optimization n>ode is considered to 

perform the recycling function; 
20 Fig. 1e is a process flow diagram complementary to the diagram of Fig. 

Id, representing tiie method steps associated witii a reclaiming function provided 

in the preferred embodiment of the optimization method according to the 

invention; 

Fig. 2 is a process flow diagram defining a sequence of method steps 
25 included in the additional pncxiuct data-generating step represented on Fig. 1c; 

Fig- 3 is a schematic representation of a cutting layout for a given piece of 
lumber tiiat has been processed according to an optimization mode adapted to 
rip-first cutting process, showing the two opposed main faces of tiie piece of 
lumber; 

30 Fig. 4 is a schematic representation of the same piece of lumber shown in 

Fig. 3 to which the method of the present invention according to an optimization 

mode adapted to crosscut-first process has been applied; 

Figs. 5a and 5b are schematic representations of pieces of lumber 

showing annangements of subdivided piece surface sections according to a first, 
35 rip-first optimization mode and a second, crosscut-first optimization mode, 

respectively; 
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Fig. 6 represents a graphical display screen generated by computer 
software, showing a table giving characteristic data related to a plurality of 
predetermined grades that can be considered in a given simulation process; 

Fig. 7 is a schematic representation of a piece of lumber characterized by 
5 the presence of several defect areas of different types as found on the two 
opposed main surfaces of the piece of lumber, 

Fig. 8 is a schematic representation of another piece of lumber 
characterized by several defect areas appearing on opposed faces thereof, in 
which a further subsection is generated considering a single-face acceptable 
10 defect criterion; 

Fig. 9 represents a main graphical display screen as generated by 
computer software implementing the method according to the invention, showing 
piece related data; 

Fig. 10 represents a graphical display screen as generated by computer 
15 software, showing a sub-menu appearing upon activation of a "File" tab; 

Fig. 11 represents a graphical display screen generated by computer 
software, showing position information related to a plurality of ripping saws that 
may be used in conjunction with the method according to the present invention, 
to periform a primary or secondary cutting pperation, depending upon the 
20 selected optimization mode; 

Fig. 12 represents a graphical display screen generated by computer 
software, showing a sub-menu appearing upon activation of an "Operations" tab; 

Fig. 13 represents a graphical display screen generated by computer 
software, showing a sub-menu appearing upon activation of a "Display" tab; 
25 Fig. 14 represents a graphical display screen generated by computer 

software, showing the list of parts to be cut according to a specific cutting bill 
entered into the computer, giving infonmation concerning each part of the list; 

Fig. 15 is a schematic representation of a piece of lumber section on 
which a selected product has been placed according to the optimization process 
30 of the invention, showing surrounding areas likely to receive additional parts 
according to a recycling approach; 

Fig. 16 is a representation of a graphical display screen as generated by 
computer software provided by the present invention, showing the result of the 
first data translating step set forth in Fig. 2; 
35 Fig. 17 is a representation of the graphical display screen as generated 

by computer software perfonning the method according to the invention, showing 
the result of the second data translating step set forth in Fig. 2; 
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Fig. 18 is a representation of the graphical display screen as generated 
by computer software perfomiing the method according to the invention, 
showing the result of alternate translation earned out at a last step set forth in 
Fig. 2; 

5 Fig- 19 represents a graphical display screen generated by computer 

software, showing the layout of selected parts to be cut according to the selected 
anangements of parts maximizing the yield value along with specific infonnatton 
related to each selected part; 

Fig. 20 represents a graphical display screen generated by computer 

10 software, giving a list of parameters related the current simulation process, 
including resulting yield data for each grade considered; and 

Fig. 21 represents a graphical display screen generated by the computer 
software, showing a table containing production data related to each part 
included In the cutting bi«. 

15 Detailed description of the preferred embodiments 

Referring now to Fig. 1a, a method of optimizing a layout of selected 
parts to be cut in accordance with the present invention will now be described in 
detail, in the context of a preferred field of application concerning lumber 
processing. However, it is to be understood that the method according to the 

20 present invention can also be advantageously employed for perfonming cutting 
optimization in other manufacturing contexts involving various types of raw 
material, such as paper, glass, fabric, plastic or rubber materials, whenever 
cutting patterns with respect to orthogonal reference axis are contemplated, 
along with geometric and grade diaracteristics of the set of parts to be cut. It is 

25 also to be understood that although the particular example to which the following 
detailed description refers involves elongated pieces of lumber such as boards, 
other types of wood raw material can be the object of the optimization method 
according to the invention, such as plywood or other kinds of wooden pressed 
panels. The two-axis optimization method according to the invention allows 

30 finding the best arrangement of parts within a layout adapted to match a specific 
piece of lumber having predetenmined known dimensions and characterized by 
the presence of defects, whose location as well as type are also known in the 
form of data stored in computer memory. The parts to be cut all originate from a 
processing bill or order, representing production requirements determined from 

35 either manufacturing of commercial needs as defined, for example, by the 
production plant manager. The main function of the optimization method is to 
identify an anrangement of parts complying with requested quantity as set out by 
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the cutting bill while seeking to recover residual raw material available between 
the areas containing the parts to be cut. The optimization method preferably 
provides relevant infomiation concerning the result of simulation processes, 
including raw material and/or economic value yield data. The optimization 
5 method according to the invention may be implemented either as part of 
production equipment such as used in rough-mills or as a stand-alone computer 
system capable of performing simulation processes to provide further flexibility of 
use. According to the former implementation, several of the optimization 
parameters can be predetemnined, therefore requiring minimal data input and 

10 control tasks by an operator. The optimization method of the invention may be 
readily embodied in the fomn of a program code as computer software adapted to 
mn on known computer hardware readily available in the marketplace, such as a 
personal computer provided with Pentium II 500 MHz or equivalent, with 128 
Mbytes of RAM using Windows NT ™ 4.0 operating software. 

15 Turning now to Fig. 1a, the optimization method according to the 

invention is capable of determining a layout of selected parts to be cut from a 
piece of raw material with respect to first and second orthogonal reference axis, 
from data representing geometric and defect-related characteristic of the piece 
as represented by data block 30 as well as from data representing geometric and 

20 grade characteristic of a set of parts to be cut represented by data block 32. 

A shown in Fig. 3. the respective outlines of first and second opposed 
main surfaces 34, 36 of a piece of lumber 38 are set at respective positions near 
or at the origin of a reference coordinates system fonmed by first and second axis 
Y and X represented at 41 and 43. By definition, the dimensions of a given piece 

25 of lumber 38 conrespond to the smallest dimensions that cover all of the parts 
forming the piece of lumber, including defect areas. It can also be seen from Fig. 
3 that each of faces 34. 36 of piece of lumber 38 is characterized by the 
presence of several defects 40, the type of which is represented by various 
graphical textures according to a predetermined graphical code tiiat will be 

30 explained later in more detail. 

Turning back to Fig- la, the optimization mefliod is initiated by a current 
optimization mode setting at an initial step 62. While the method can be 
implemented considering a single optimization mode such as basic rip-first or 
crosscut-first mode, the computer software Is preferably programmed to allow the 

35 selection of a desired optimization mode among a plurality of optimization 
modes, including rip-first crosscut-first, single-length/all-width crosscut-first 
(panel cutting), alMengtii/single-width rip-first, as well as any combination of 
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these basic optimization modes such as hybrid and comparative yield 
optimization modes. As shown in tiie example of Fig. 3, tiie rip-first optimization 
mode considers first dimension values defined with respect to a first axis Y at 
41, according to which respective width of series of aligned parts 42 and 44 are 
5 defined. For example, parts 42, 44 can be flooring wood products that are 
characterized as having a common, standard width (ex. 3 inch) without being 
associated with specific length requirements from the cutting bill, so that ttie 
computer software may determine a particular value for each wooden member 
42, 44 in order to maximize either raw material or economic value yield. It can be 

10 appreciated from Figs. 3 and 5a. accorting to a typical rip-first optimization 
mode, tiiat the elongate piece of lumber 38 shown has its longitudinal axis 46 in 
parallel to the second reference axis X , while the same piece of lumber 38 has 
its transverse axis 48 parallel to first reference axis 7 . 

In an analog manner, for a typical crosscut-first process, it can be seen 

15 from Fig. 4 that each panel included in a given subsection containing series of 
pane\s 54, 55, 56, 57 and 58 has a predetermined standard lengtii. while the 
width dimension of each panel does not correspond to any requirement from the 
cutting bill, tiie width value for each panel being determined according with 
highest yield criteria as will be later explained with more detail. In the aosscut- 

20 first example of Fig. 4, it can be seen tiiat tiie elongated piece of lumber 38 has 
now its longitudinal axis 46 parallel to a first reference axis Y at 41', while having 
its transverse axis 48 parallel to a second reference axis Xai 43'. The rip-first 
cutting process is furUier characterized in that the primary cuts represented by 
cutting axis 50 in Fig. 5a extend parallel to the longitudinal axis 46 of the piece of 

25 lumber 38. As will be later explained in more detail, the secondary cutting 
operation in a direction parallef to transverse axis 48 of the elongated piece of 
lumber 38 shown in Fig. 3 is performed sequentially to the primary cut operation 
along axis 50 shown in Fig. 5a in order to produce series of selected parts 42, 44 
illustrated in Fig. 3. The crosscut-first mode of optimization involves a cutting 

30 process approach according to which the first dimension values considered in 
priority are defined witii respect to a first axis Y as shown at 41' in Fig, 4, 
wherein the elongated piece of lumber 38 has its longitudinal axis 46 parallel to 
first reference axis Y at 41', while having its transverse axis 48 parallel to 
second reference axis X at 43'. The crosscut-first cutting process is furttier 

35 characterized in that the primary cuts represented by reference axis 52 in Fig. 5b 
are performed first in a direction parallel to transverse axis 48 defined by piece of 
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lumber 38, followed by secondary cutting operations in a direction parallel to 
longitudinal axis 46 defined by piece of lumber 38 shown in Fig- 4 to produce 
series of parts 54, 55, 56, 57 and 58 as wilt be explained later in more detail. 
Figs. 3 and 4 illustrate general cases involving all-length, all-width rip-first and 
5 crosscut-first optimization modes, respectively. The panel cutting optimization 
mode is a species of crosscut-first mode explained above, in which the parts to 
be selected are characterized by having a single, common lengtii while having 
various widths, so that the primary cutting operation is to be perfomned over the 
whole widtii of a given piece of lumber. The alHength/single-widtti cutting 

10 process is a particular case of the basic rip-first process explained above, in 
which the selected product are characterized as having a single, common widtii 
while having various lengths, and the primary cutting operation is to be 
performed over the whole length of a given piece of lumber. Hybrid modes of 
optimization involve more than one basic mode of optimization to detemnine the 

15 anrangement of parts maximizing the yield value associated wiUi the optimized 
layout to be produced. For example, such hybrid optimization mode may 
generate an optimized layout consisting of an arrangement of parts such as 
panels 42, 44 as well as components 60 having different dimensions as shown in 
dotted lines in Fig. 3. As opposed to panels 42, 44 shown in Fig. 3 and flooring 

20 wood members 54, 55, 56. 57 and 58 shown in Fig. 4, each component 60 
shown in Fig. 3 has dimensions in tiie directions parallel to axis Y and X 
corresponding to a requirement of the cutting bill according to which the 
optimized layout is to be generated. Preferably, for hybrid optimization modes, an 
optimization time limit for processing each piece will be established and shared 

25 between all optimization modes Involved. 

Tuming back to Fig. la, the current optimization mode setting step 62 
convenientiy initializes an optimization process by defining a reference 
coordinate system adapted to the selected optimization mode on the basis of a 
reference axis assignment prindple according to which a)ds Y and X are 

30 arbitrarily chosen to be first and second ortiiogonal reference axis to perform rip- 
first optimization as applied in the example of Fig. 3, while same axis X and 
Fare chosen as first and second orthogonal reference axis for the purpose of 
crosscut-first optimization as applied in the example of Fig. 4. Such pemnutation 
of reference axis allows tiie use of essentially tiie same algoritiim to perform 

35 optimization processes of the invention according to two or more optimization 
modes. 



wo 2004/051523 



PCT/CA2003/001872 



12 

Turning again, to Fig. 1a, the method further comprises a step 63 of 
defining a subset of the set of parts listed in the cutting bill, characterized by a 
predetemiined grade value / and associated with a predetermined group of first 
dimension values defined with respect to the first axis, namely with values 
5 defined with respect to axis Y in the case of a rip-first optimization process, as 
illustrated in Fig. 3. In the application of the method of the invention for pieces of 
lumber, the value assigned to grade i specifies a predetermined level of wood 
quality, each grade i being associated with one or more criteria related to 
characteristics of the raw wood material that will be used to produce a given 

10 product to which that specific grade i has been assigned. Preferably, for a 
number n of predetermined grade values that can be assigned to each product 
to be cut, an increasing cutting priority order is considered to allow 

systematic generation of the an^ngement of parts offering the best raw material 
or economic value yield. According to the invention, the optimal cutting layout 

15 solution is sought on the basis of subsets that are defined from a set of parts to 
be cut as established by the cutting bill defining the optimization mode and grade 
currently processed. In the example shown on the display screen of Fig. 6, the 
table 64 shown presents the criteria characterizing six (6) distinct grades of 
decreasing quality listed in optimization priority order, namely clear, clear single- 

20 face, superior, dark, thin and reclaimed, each of which grades being associated 
with a distinct combination of defect features related to the presence or absence 
of knot, wane, coloration and missing wood. The given priority order assigned to 
each grade detennines the placement sequence of corresponding parts. For 
example, while the presence of any of the above-mentioned defects on the raw, 

25 wood material section would not be usable for a product to which a dear grade 
has been assigned, a raw wood material section characterized by only coloration 
without the presence of other kind of defects will always satisfy the requirements 
of either dark, thin or reclaimed grades, and may satisfy the requirements of 
either clear single-face or superior grades provided the coloration feature is 

30 found on a single face of the considered section of raw wood material. In order 
to allow accurate primary and secondary cutting operations while having the 
option of selecting parts of higher grade, the algorithm used by ttie computer 
software to perfomn the optimization metiiod of the invention preferably considers 
both faces of a processed piece of lumber in the determination of whether a 

35 single-face defect should be considered as acceptable or unacceptable 
depending upon the specific grade considered, as will be explained later in more 
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detail with reference to Figs. 7 and 8. It can be seen from table 64 of Fig, 6 that 
to each grade can be associated a specific colour/shade code that could be 
conveniently used to facilitate the identification of parts listed in the given cutting 
bill that are assigned to a same grade, as will be explained later in more detail 
5 with respect to Fig. 14. 

Turning now to Fig. 9 showing a main display screen generated by the 
computer software, relevant data representing geometric and defect-related 
characteristics of a given piece of lumber such as board no. 6 In the instant 
example is presented graphically. The main display screen of Fig. 9 provides 

10 various types of information, namely identification number assigned to the board 
for which data is displayed; board dimension data including length, width and 
surface area; initial and final clear-out positions relative to the origin of the 
reference coordinate system used, which values are chosen to ensure that the 
dear-out end portion of the board mainly includes unusable raw wood material; 

15 and number of identified defects found on either main faces of the current board. 
Furthenmore, the main display screen preferably includes a table shown at 59 
giving data specific to each listed defect, such as assigned identification defect 
number, position coordinates of each defect expressed as a pair of coordinates 

(Xmin ' y^^' ^Xmax' ^ ^Nhlch roprosent the smallest rectangular area that 

20 completely contarn the defect; identification of the face on which the defect 
appears; arid type identification of the defect, namely knot, wane, coloration or 
missing wood as indicated in a legend generally designated at 66 on Fig. 9. The 
computer software implementing the method is configured so as tiie main display 
screen of Fig. 9 furtiier includes a schematic plan view of both first and second 

25 faces 34, 36 of the current board into which tiie listed defects are represented by 
corresponding rectangular areas using the appropriate graphical code given by a 
legend 66. In the example shown, a small knot 68, larger knot 70 as well as 
coloration areas 72, 74 are associated with first face 34, while knot 76 and 
coloration areas 78, 80 are associated with second face 36. Finally, at the 

30 bottom portion of the display saeen of Fig. 9, appropriate scales for axis Y and 
X used to represent first and second board faces 34, 36 are given. 
Convenientiy, data regarding the definition of all assignable types of defects is 
contained in a configuration file to which the computer software has access. The 
types of defects that could be considered can be defined by anyone skilled in the 

35 art according to the specific application involved. In the context of lumber 
optimization application, typical defects include: edge-related defects associated 
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with missing wood at an edge portion of a piece of lunnber generally due to edge 
defomnation or natural bending of the piece; wane characterized by missing wood 
generally located at edges of a piece of lumber mainly due to the generally 
circular cross-section of timber or to over-cutting operation at the planner stage; 
5 missing wood at the surface of a piece of lumber generally due to over-cutting 
either at debari<ing, pre-cutting or planning stages; and unacceptable dark knot 
and coloration areas that are unusable for specific application such as furniture 
manufacturing. 

Tunning now to Fjg. 10, the pull-down sub-menu 82 upon activation of the 
10 "File" tab 84 gives access for the operator to an "Open optimization process" 
function as indicated at 86 that allows to open a simulation data file that has been 
previously stored in memory, the activation of which function will clear any 
cunnently downloaded data to substitute therefor data related to the selected 
simulation file. Preferably, each simulation data file refers to a number of 
15 subordinate data files having a different identification extension, such as lists of 
raw material pieces data file, cutting bill data file, grade definition data file, 
simulation result data file and replacement parts data file, the latter type of data 
being used to replace a part having been selected in a layout, by a new required 
part listed in a cunrent cutting bill. A typical simulation data file includes data 
20 related to cutting bills, a list a raw material such as pieces of lumber, all 
simulation parameters as well as information regarding selected grades. Table 1 
gives a list of typical simulation parameters with exemplary values contained in a 
simulation data file. 
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QIIUIIII ATinM DAPAMPTPP ' 


\/AI IIP 


iName or simLfiauon 


o im u na uon- 1 


Simulation start index (pieces data file) 


1 


oirnujation ena inaex (pieces aaia vie) 


c 
O 


oimuiauon moae 


Kip-Tirsr 


Secondary simulation mode 


none 


Tertiary simulation mode 


None 


Minimum yield reference value (%) 


50 


Units system 


Metric (mm) 


Piece thickness 


on 
20 


Optimization time limit (s) 


600 


Weighing type 


Reclaimed and factor 


Weighing 


Yes 


Weighing factor- panel 


0.75 


Weighing factor - flooring wood 


0.75 


Minimum width — panel 


13 


Maximum width — panel 


175 


Minimum length — flooring wood 


254 


Maximum length — flooring wood 


21 o4 


Reclaiming function 


Yes 


Minimum reclaimed length 


4oo 


Maximum reclaimed length 


HOOT 


Minimum reclaimed width 


12 


Maximum reclaimed width 


4o2 


Ripping system length 


50o 


Lrrosscuuing system lengtn 


oU4o 


A||_lannfh vali i^q 

AMl~lCl|^UI VcllUCo 


1 C?o 


All-width values 


Yes 


Maximum remaining length 


25 


Saw kerf on edges 


Yes 


Ripping saw kerf width 


4 


CrosscutHng saw kerf width 


4 


Replacement rule 


Rule 1 
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AII-prioiKy allowed 


Yes 


Replacement length threshold " 


51 


Grade optimization 


Yes 



TABLE 1 



Conveniently, if the current downloaded simulation file is afterward modified, 
the computer software will cause the display of a store confinr^ation dialog box to 
allow the operator to store updated data before leaving the session. Convenient 
5 functions such as "Store". "Import file", "Export file". "Print file", "Print whole 
optimization", "Preview", "Printer configuration" and "Quit" accessible through sub- 
menu 82 are preferably provided, such basic functions being readily implemented 
into a computer software through conventional programming. In particular, the "Prinf 
may be used to print a cunrently displayed screen showing data regarding a cunrent 

10 piece of lumber, while the "print whole optimization" funcBon may be used to print all 
optimization layouts downloaded in memory associated with all pieces of lumber 
considered by the optimization process. 

Refem'ng now to Fig. 11, there is represented a graphical display screen 
generated by the computer software, presenting posftion infomnation related to a 

15 plurality of ripping saws that may be used to perfonn either the primary cutting 
operation in the context of a selected rip-first optimization mode, or a secondary 
cutting operation when a crosscut-first optimization mode has been defined. Sawing 
system characteristics such as length and number of saws are shown, as well as a 
table 105 giving infomnation related to each saw, namely saw identification number, 

20 position of saw (origin), width of cut, saw kerf width and indication whether each 
saw is mobile or fixed. A graphical representation of the position of each saw is 
presented within a firame 107. 

Referring now to Rg. 12, a furtiier pulkJown sub-menu 88 generated upon 
activation of an "Operations" tab 103 is shown, which sub-menu 88 gives access for 

25 tiie operator to a group of ftjrttier functions including a first "Start simulation" function 
at 90 allowing the operator to start a simulation process using the currently 
downloaded data contained in tiie selected simulation data file. A "New simulation 
process" function designated at 92 allows to clear the data regarding Oie cunient 
simulation process as vyell as ttie simulated part data defined by the current cutting 

30 bill to display a new blank screen, allowing the operator to initialize a new simulation 
process. The "Simulation series" function designated at 94 allows the operator to 
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cause the computer software to perform without interruption a predetemiined series 
of simulation processes between which specific parameters values are downloaded 
according to a pre-programmed control "batch" file, referring to sub-files such as 
simulation data files, simulation result data files, and production data file when the 
5 cutting bill is associated to a specific plant production. Each listed file is downloaded 
before the execution of each simulation process, and a conresponding series of 
simulation result files are generated and stored in memory. 

Tuming to Fig. 13, a further pull-down sub-menu 109 generated upon 
activation of a "Display^ tab 111 is shown, which sub-menu 109 gives direct access 

10 to functions such as tool bar and state bar setting functions, main board screen as 
shown in Fig, 9, as well as to otiier displayed screens that will be described later in 
detail, namely, simulation parameter, grades, cutting solution, results and enx)r 
screens. Conveniently a Tools" tab at 113 may be provided to give access to any 
optional data analysis functions that could be implemented into the computer 

1 5 software if needed. 

Referring to Fig- 14, the cutting bill display screen gives the number of length 
values involved by the current set of parts to be cut whose characteristics are 
presented in the fomn of a table 115, containing the following data for each listed 
part length value, width value, product type, requested quantity, remaining quantity 

20 (units), remaining quantity in %, grade, production identification number, part code, 
replacement priority indicator, weighing factor value, price (for economic value- 
based yield), active set of products indicator and infinite quantity indicator, the latter 
parameter being used whenever a limit requested quantity is not assigned to a 
particular part to be produced. Conveniently, the data related to each part of a 

25 cunBntly active set of parts is presented in dark characters, while data related to 
remaining, replacement parts is presented in grey shade. The part replacement 
function and assodated rules \Anll be explained later in more detail. The computer 
software preferably provides other functions, such as to allow part addition into a 
cunnent cutting bill, to allow parameter editing on several parts simultaneously, to add 

30 grades for the parts, to sort data in tiie various display screens or to print data 
reports. 

Tuming back to Fig. la, a group of first dimension values as set at step 63 
allows the determination of simulation parameters characterizing each simulation 
trial as will be now explained on the basis of the following example, wherein a 
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wooden board having a width of 150 mm is subjected to a rip-first optimization 
process considering a subset of three (3) distinct width values, namely 100 mm, 50 
mm and 25 mm. For each trial, a numerical vector is generated as represented in 
Table 2, wherein to each column is associated a specific index designating one of 
5 the first dimension values of the group established from the subset of parts from 
which an arrangement is defined at each trial. 



Indexl 


Index2 


IndexS 


Index4 


IndexS 


Index6 


1 


2 


-1 


-1 






1 


3 


3 


-1 






2 • 


1 


-1 


-1 






2 


2 


2 


-1 






2 


2 


3 


3 






2 


3 


3 


3 


3 




3 


1 


3 


-1 






3 


2 


2 


3 






3 


2 


3 


3 


3 




3 


3 


3 


3 


3 


3 


TABL 


E2 



The dimension of vector is detemriined according to the following relation: 



VectorDim^PieceWidtl^^^^^^^^^ (1) 

15 

wherein: 

PieceWidth is the width dimension of the piece; 

SmallWidth is the smaller width dimension of the group of first dimensions; 

and 

20 KerfWidth is the width of the kerf produced by the ripping saw to be used. 

The VectorDim represents the maximal number of piece sulxlivisions that 
can be obtained for the optimization trial solution, through a primary rip-flrst cutting 
operation with respect to second reference axis X represented on Figs. 3 and 5a, 
considering the selected optimization mode and grade. In the example given in 

25 Table 2, VectorDim is equal to six (6), i.e. fndexl to Index6, and ten (10) different 
optimization trial solutions have been generated according to a matrix filling 
algorithm that will be now explained. At the beginning of each iteration step, all 
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matrix elements of the current matrix are filled with an initial value that has been 
arbitrarily chosen as "-V to indicate that no conresponding number of subdivisions, 
characterized by the width dimension associated with the corresponding index, has 
been initially assigned to the index parameter. The index values are obtained 
5 through the first iteration step according to the following operation sequence: 

CurrentWidih=Piecemdth 
FROMi=0 to i^NumWidth 
{ 

1 0 Width^GetWidthfiJ-^Kerfrndth 
IF Width<CurrentWidth 

{ 

MatrixIndex=0 
15 FOR each R^eat 

{ 

CurrmtWidth^CurrentWidth'Width 
MatrixjMatrixIndex] =1 
MoOrixIndex-MatrixIndex'i'l 

20 ; 

; 

; 

wherein: 

25 CurrentWidth is a parameter representing the current value of the remaining 

piece width; 

NumWidth is a parameter representing the number of width dimensions 
contained in the group of first dimension values; 

Width is a parameter representing the effective width dimension value 
30 corresponding to a current index of the matrix; 

GetWidth(i) is a parameter representing the width value conresponding to 
width dimension of index /; 

Repeat represents the number of subdivisions associated with a cunrent 
index of the matrix; and 
35 Matrixlndex is a parameter designating a current index. 

Subsequent matrix line filling steps are performed by finding a last index that 
can be decomposed, and by substituting for such last index a following index value 
corresponding to the next smaller width dimension that will t>e used to fill the 
remaining available piece surface. It can be appreciated that index values z=l,6 
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shown in the exannple of Table 2, are assodated with width dimensions of 
decreasing values. Such last index finding step can be performed according to the 
following position finding sequence: 



5 DecompIndex=0 

FROMi=^0 to i==^VectorDim 

{ 

IFMatrix[IJ<SmalimdthIndex 

{ 

10 DecompIndex=I 
} 

} 



wherein: 

15 Decomplndex is a parameter designating an index value corresponding to 

the current position to be decomposed; and 

SmallWidthlndex is a parameter designating an index value corresponding to 
the smallest width of the group. 

Then, the matrix must be further modified in order to decompose the position 
20 found, by first recalculating the remaining available width until the position just 
preceding the matrix position to be decomposed is reached. Then, the selected 
position is decomposed while the remaining cells of the matrix are filled by 
calculating width of used material, according to the following sequence: 



25 CurrentWidth^PieceWidth 

FROMi=0 to i==(Deco7npIndex-I) 
{ 

CtirrentWMth=Currentmdih'GetWidth(MatrixflJ) 

} 

30 

Index-Matrixp)ecompIndex 
Index^Index+1 
Matrix[DecompIndex] -Index 

Cunentmdth=CurrmtWidth'GetWidth(MatrixlpecompIndex]) 

35 

FROMj=0 to j^NumWidth 
{ 

mdth-=GetmdthO)+KerfPFidth 
IF Width<CwTentWidth 
40 { 

R^=CurrentWidt1^^^^ 
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MatrixIndex=DeconipIndex+l 

FOR each Repeat 

{ 

CurrentWidth^CurrentWidth'Wtdth 
5 Matrix[MatrixIndexJ=I 

MatrixIndex=MatrixIndex+l 

} 

} 

} 

10 FROMk==0 to h=MatrixDimension 
{ 

Index=Matrix[k] 

IF IndexOSmallWtdthlndex 

{ 

15 FilledMathx^ FALSE 

} 

} 

FiUedMatrix= TRUE 



20 wherein: 

Index is an incremental parameter, 

MatrixDimension designates the specific dimension of the matrix which is 
equal to the resulting value of nxn ; and 

FilledMatrix is a flag indicating whether matrix filling is completed or not, 

25 The numt>er of different lines fomning a current matrix defining the 

optimization trial solutions to be generated is dictated by the PieceWidth value 
characterizing the currently processed piece, the cutting layout of which has to be 
optimized, sis well as by the number of width dimensions fonming the group 
characterizing the subset of parts serving as a basis for the optimization trial The 

30 optimization process is completed when all trials have been perfonmed so that all 
cells of the last trial correspond to the index value of the smallest width as shown in 
Table 2. It is to be understood that the foregoing algorithm for building the trial 
matrix as applied to the rip-first optimization mode can be readily adapted to other 
types of optimization modes such as crosscut-first or hybrid optimization modes by 

35 providing appropriate reference axis system permutation as explained before 
regarding the optimization mode setting step 62 shown in the flow diagram of Fig- 
la, For example, for an application involving crosscut-first optimization mode, the 
length dimension of the pieces would be substituted for the vwdth dimension referred 
to in the foregoing example involving rip-first optimization mode, while the 

40 subdivisions of each piece would be established from a group of length values as 
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first dimension values defined with respect to axis Y considered as the first axis as 
shown In Figs. 4 and 5b. It is pointed out that the foregoing algorithm allows to 
generate data defining an arrangement of subdivided piece surface sections to be 
obtained through a primary cut operation, corresponding to step 65 in Fig. la, which 
5 piece surface sections are designated at numeral 67 and 69 in Figs. 5a and 5b 
respectively associated with rip-first and crosscut-first optimization modes as 
described above. It is also pointed out that at this stage of the method, none of the 
optimization trial solutions includes any subdivided piece surface section having 
width/length dimensions characterizing parts to which are assigned any grade that is 

10 different from the selected grade for the current optimization process. At a later 
stage of the optimization process, the l>est optimization trial solution will consider 
other grades of lower quality according to a recycling approach that will be later 
explained in detail. The data generating step 65 is performed from the matrix as 
filled with width indicative index values taken from the group of width values 

15 considered as first dimension values, wherein to each width value corresponds a 
number of parts characterized by that common width value as defined in the current 
cutting bill. 

Turning back to Fig. la, following step 65 is a further step 71 of generating 
data finding one or more subsections included in each piece surface section 

20 generated at preceding step 65, according to the geometric and defect-related 
characteristics of the piece of luml)er under processing. Each of the optimization trial 
generates a number of piece surface sections to be obtained through a primary 
cutting operation, each of which piece surface sections being afterwards sut>divided 
independently to produce a basic cutting optimization solution as will be later 

25 explained in more detail. Such basic cutting of the piece surface sections on the 
basis of the data generated at step 71 preferably considers a two-face lateral 
exploration of each piece of lumber, wherein all defect identified is sorted so as to 

obtain a defect list in temns of coordinates r ■ and v - 

Tuming now to Fig. 7, a typical piece of lumber section is schematically 
30 represented by top view of first face 73 and bottom view of second face 73'. It can 
be seen that respective positions of secondary cutting lines 75, 77 and 79 are 
detemnined so that each of the three (3) subsections generated by these secondary 
cutting lines, as designated at 81, 83. 85 on first face 73 and at 81', 83', 85' on 
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10 



15 



second face 73', is either free of any defect or contains one or more acceptable 
defects depending on the current grade, as explained above with reference to Fig. 6. 
Furthermore, each subsection has a length larger than the minimum length 
characterizing the parts contained in the considered subset and has a width value 
that is equal to the width of the whole section whose faces 73, 73' are represented 
on Fig. 7. The data-generating step 71, along with series of steps 87. 89, 91, 93 
and 95 which are necessary to optimize the cunrent subsection according to the set 
optimization mode as explained later in more detail, can be implemented into the 
computer software according to the following algorithm: 

NumDefect=GetNumDrfect(Piece) 

XMax^XSectionMax 
XMin^XSectionMinSelectFace^l 

FROMi=^0 to i='NumDefect 
{ 

CurrentDefect=^GetDrfect(i) 
20 IF CurrentDefect intersects the current piece section 

{ 

Def€ctCkttegory~GetDefectCategory(CiurentDefect,SelectGrade, 
SdectFace) 

25 IFDefeciCategoryoAcceptable 
{ 

AycdlLength=XMinIntersection-Xmin 
IF AvailLengh>MmLengfh 

f 

30 Generate the Subsection 

Optimize the Subsection according to set optimization mode 

} 

IF SinsleFaceDefec~mUE 
{ 

35 XMin^XMinD^ect 

IF SdectFace==l 
{ 

SelectFace=2 

} 

40 ; 

IFNOT 

XMn=XMaxDefect 

} 

} 

45 ; 
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AvcalableLength=XMax-XMin 

IF AvailableL€ngth>MiiiLength 
5 { 

Generate the Subsection 

Optimize the Subsection according to set optimization mode 

} 

10 wherein: 

NumDefect represents the number of defects identified for a currently 
processed piece of lumber, 

XMax designates the maximum coordinate value along second axis X in 
relation with a currently processed subsection; 
15 XSectionMax designates the maximum coordinate value along second axis 

X in respect of a currently processed section; 

XMin designates the minimum coordinate value along second axis X in 
relation with a currently processed subsection; 

XSecSonMin designates the minimum coordinate value along second axis X 
20 in respect of a currently processed section; 

SelectFace designates a selected one of first and second faces of the piece 
of lumber; 

CurrentDefect designates the defect cunrently processed; 

GetDefect(i) designates defect-characterizing data associated with a defect 

25 /; 

DefecfCafego/y designates a currently processed defect category; 

GetDefectCategory designates category-characterizing data according to the 
cun^ent defect, selected grade and selected face; 

SeledGrade designates the selected grade under processing; 
30 AvailabfeLength designates available length value of the currently processed 

subsection; 

XMinlntersection designates the coordinate value along second axis 
X where the cunrent defect intersects the currently processed section; 

MinLength designates the minimum length value characterizing the subset of 
35 parts considered; 
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SingleFaceDefect is a logical parameter indicating Nvhether a single-face 
defect elimination is required; 

XMinDefect designates the minimum coordinate value along second axis 
X related to current defect under processing; and 
5 XMaxDefect designates a maximum coordinate value along second axis X 

in respect of a cunrent defect under processing. 

Turning again to Fig. la, data generating step 87 by which, for each 
subsection, a plurality of arrangements of parts selected from the subset of parts to 
be obtained through the secondary cutting operation, will be now explained in 

1 0 greater detail. The ultimate purpose of data generating step 87 is to fill a defect-free 
subsection as generated at prior step 71 , by placing therein the arrangement of parts 
selected from the previously established subset of parts that will maximize surface 
utilization of the cunrently processed section. For so doing, a plurality of distinct 
arrangement of parts must be considered and qualified, using a matrbc-based 

1 5 approach that is similar to the one applied for generating optimization trials obtained 
at step 65. A trial matrix is generated through position penmutation of selected parts. 
It is pointed out that as opposed to the trial generation algorithm explained above in 
respect of step 65, placement order within a same subsection does not have any 
effect on the resulting raw material yield. Preferably, each arrangement of parts is 

20 generated so as to maintain at least one product of each width value that has been 
already selected, to preserve parts having the largest width that have been selected 
first. For example, with a subsection of ^00mm width and 400mm length, lets 
assume that the following subset of parts has been established according to table 3: 



Part 


Width 


Length 


Number 


mm 


mm 


1 


120 


350 


2 


100 


200 


3 . 


90 


150 


4 


75 


150 


5 


50 


100 



25 TABLE 3 
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Each line of the secondary optimization trial matrix shown in Table 4 can be 
filled with appropriate index values, wherein only two lines are generated in the 
instant example considering the above-mentioned condition according to which 
decomposition is canied out while ensuring that at least one product of each width 
5 that has been already selected is included within the annangement. 



Indexl 


Index2 


Index3 


Index4 


2 


2 


-1 


-1 


2 


3 


-1 


-1 


TAB 


LE4 



The data generating step 87 may be perfonned using the following algorithm: 

10 

WHILE LastTrialCompleted=FALSE 
{ 

Siaface^O 

LastTrialCompleted^FIND.DECOhfPINDEX 

15 

FnX.MATKIXSECOmCUT 
NtmiPart=GetNumPqrt(Subsection) 

FOR i=0 to i^NumPart 
20 { 

Part^GetPart(MatnxIndex[lJ) 
Swface^nSwface-^GetSiirfaceCPart) 

} 

} 

25 

wherein: 

LastTrialCompleted is a flag indicating whether a last simulation trial has 
been completed or not; 

Surface is a parameter whose value represents the surface area obtained 
30 from the current annangement of parts associated with a currently processed 
subsection; 

FIND.DECOMPINDEX designates a subroutine the purpose of which 
consists of finding the value of a cunrently decomposable index and returning a flag 
value as an indication of the end of the decomposition process; 
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FILLMATRIXSECONDCUT designates a subroutine the purpose of which is 
to fill the optimization trial matrix associated with the currently processed subsection; 

NumPart is a parameter whose value represents the number of parts 
contained in the arrangement associated with a current optimization trial; 
5 GetNumPart(Subsection) designates a function according to which the 

number of parts included in a cunrently generated arrangement is assigned to the 
NumPart parameter; 

Part is a parameter providing identification of a specific part contained in the 
arrangement of parts associated with a current optimization trial; 
10 GetPart(Matrixlr}dex[i]) designates a function according to which a value for 

the part parameter is found corresponding to the ojrrently processed Matrixindex P] 
value; and 

GetSurface(Part) designates a function according to which the surface value 
associated \Anth a specific part is found for the purpose of surface calculation. 
1 5 . The FIND.DECOMPINDEX subroutine may be implemented according to the 
following algorithm: 



IFMatrixfOJ=--J 

Return FALSE 

20 

FROMi=-0 to i=VectorDim 
{ 

Ifidex=GetIndex(i) 
IF IndexoSmaJlWidtMndex 
25 { 

IF LastPiece=FALSE 
{ 

Deco7npIndex=i 
Return FALSE 

30 ; 

; 

; 

wherein: 

35 Getlndex(i) designates a function according to which the index value 

corresponding to the current value for / is found and assigned to the index 
parameter. 

It can be appreciated that according to the at)ove algorithm, whenever the 
matrix value at first position "0" responds to a "0" value, it immediately follows that 
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the current matrix has not been previously initialized, and the return flag value 
"FALSE* causes the algorithm to generate a first line within the current matrix. The 
algorithm will retum a "TRUE' flag value whenever there is no further matrix position 
index to be decomposed and will retum "false" flag value whenever there is one or 
5 more further decomposable index positions. Whenever the condition Matrix[0]= -f is 
satisfied, the first matrix line can be generated as follows: 

Ct(rrentLengfh=SubsectionLength 
Position=0 

FROMi^PartStarandex to i^PartEndlndex 
{ 

Part=GetPart(i) 

PartLength=GetLength(Part)+ KerfWidth 
PartWidth^GetWidth(Part) 

WHILE (PartLength<CurrentLength) AND (P'artWidlh<SubsecnonWidth) 
{ 

CurrentLength=CurrentLength-PartLength 
Matnx[Position]-i 
Position=Position+] 

} 

} 

wherein: 

. 25 CurrentLength is a parameter representing the current value of the available 

subsection length; 

SubsecUonLength represents ttie total length value of the cunnent subsection; 
Subsecb'onWidth represents the total width value of the cunnent subsection; 
PartStartlndex represents a start index value associated with the subset of 
30 part^ selected according to the current grade and optimization mode; 

PartEndlndex represents an end index value associated with the subset of 
parts selected according to the current grade and optimization mode; 

GetPart(i) represents a function according to which data corresponding to a 
part / selected from the subset of parts is found; 
35 PartLength represents an estimated value for the resulting length of a part / 

considering the kerf width value; 

PartWidth represents the width value associated with part / ; 
GetPartWidth represents a function according to which the current width 
value of part / is assigned to PartWidth parameter. 



10 



15 



20 
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40 



Once the first line of the subsection trial matrix has been generated, one or 
more further lines of index values can be generated as part of the 
FILLMATRIXSECONDCUT subroutine according to the following algorithm: 



CtarentLength=^SubsectionLength 



FROMi=^0 to i=(DecompInd€x-l) 
{ 

10 ParP=GetPart(i) 

PartLength^GetLengfk(Part)+ KerfWidih 

CiirrentLmgfh=CtinentLmgth-FartLength 

} 

1 5 Index=Matrix[DecompIndex] 
Index-Index-^ I 

SelectPart=FALSE 

20 WHILE (SelectPart-^FALSE AND Index<PartEndIndex) 
{ 

Part^GetPart(Index) 
PartLength^GetLength(Part)+Kerff¥idih 

25 IF Length<CurrentLmgth 

{ 

SelectPof^TRVE 

CiirrentLmgth=CurrentLength'PartLength 
Matnx[DecompIndex]— Index ' 

30 ; 

IFNOT 
{ 

Matrix[Dec<mpIndexJ- -I 

} 

35 ; 

FROMj=DecompIndex to j-VectorLength 
{ 

Matrix[j]= -1 

} 



NextIndex=DecompIndex-^I 



FROMk=Index to k=PartEndIndex 
{ 

45 Poj^GetPartfk) 

PartLength==GetLengfh(Part)'\^KerfLength 
PartWidth=GetWidth(Part) 
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WHILE (PartLength<Ci4rrentLenph) AND (PartWidth<Subsectionmdth) 
{ 

CurrenlLength-CurrentLength-PartLength 
Matrix[NextIndex]-k 
5 NextIndex=NexiIndex-^l 
} 

} 

wherein: 

10 SelectPart designates a flag parameter Indicating whether a further part has 

been selected or not; and 

Nextlndex is an incremental parameter used by the algorithm to complete 
matrix filling for the remaining index positions of further lines. 

In the example of Table 4, the matrix shown includes two (2) lines of index 

15 values connesponding to two (2) optimization trial solutions for a current subsection. 
After having generated data defining for a current subsection, the plurality of 
anrangement of parts to be Included therein and selected from the subset of parts, 
i.e. parts 1-5 in the above example, which parts will be produced through a 
secondary cutting operation with respect to a first reference axis 7 shovkffi in Figs, 3 

20 and 5a, the method according to the invention further includes a step 89 as shown in 
Fig. 1a according to which yield values associated with all arrangements of parts are 
estimated. It can be seen from Table 4 that the first trial solution Involves the 
Inclusion of two (2) parts of index "2" assigned to Indexl and Index2, resulting to 
100% surface coverage of the cun^nt subsection (40 000 m/n^, while the second 

25 optimization trial solution for the same subsection involves the inclusion of one (1) 
part of index "2" along with one (1) part of Index "3". assigned to Indexl and Index2 
respectively, the latter arrangement of parts leaving a remaining surface area 
unused (500 mm^. The method further proceeds with a step 91 of comparing the 
yield values to select the anangement of parts having a highest yield value. In tiie 

30 example for which optimization results are presented in Table 4, the first 
optimization trial solution for the current subsection whose index values are given in 
the first line of the matrix is found to hiave tiie highest surface yield value compared 
with the yield value provided by the optimization trial solution whose index values 
are given at the second line of the matrix. Lefs now consider another example of 

35 the proposed method involving a subset of parts whose characteristics are given in 
Table 5 with the resulting subsection optimization matrix shown in Table 6. 
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Part 


Width 


Length 


number 


(mm) 


(mm) 


1 


120 


350 


2 


100 


200 


3 


100 


120 


.4 


90 


170 


5 


50 


100 



TABLES 



Indexl 


Index2 


Index3 


Index4 


2 


2 


-1 


-1 


2 


3 


3 


-1 


2 


3 


4 


-1 



5 TABLES 

It can be appreciated that the third optimization trial solution for the current 
subsection is found to be the best solution following comparison of its resulting 
surface yield (49 000 mm^) with the resulting yield value obtained with the first trial 

10 solution (40 000 mn?^) and second trial solution (44 000 mm^). Even if part no. 4 is 
characterized by a width value lower than the corresponding width value of part no. 
3, the anangement part lengths provided by the third trial solution allows complete 
filling of the available surface area of the current subserfion. It is to be understood 
that although the maximizing surface criterion has been preferably applied in the 

15 above examples, an economic-based criterion or combined economic/surface yield 
criterion may be also applied using the same basic algorithms as described above. 
Although the above algorithms are adapted to rif>-first optimization mode, these 
algorithms can be readily adapted to other optimization modes by any person skilled 
in the art of computer programming, using the axis permutation principle as set forUi 

20 before, or any other equivalent approach. For example, the axis pemnutation may be 
obtained using to the following operation sequence : 

Store=XMin 
Xmin-YMin 
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YMin^Store 
Store^YMax 
YMax^Xmax 
JOnax=Store 

5 wherein Store is a transfer variable. 



In the particular case of panel prcxluction, crosscut-first optinnization mode, the 
following algorithm may be used to perfomied step 87 of Fig. 1a for a current 
subsection: 

10 

CurrentLength =SubsectionLength 

PositionX=^XminSubsection 

AvailWidth=SubsectionWidth 

1 5 FROM i=PartStartIndex to i^PartEndlndex 
{ 

Part=GetPart(i) 
Lengih=GetLeng^h(Part) 
IFLength<CurrentLength 
20 { 

K^eat /MaxBoardWidth-^KirfWidth 

WHILE Repeat>0 

{ 

Select a Panel having width=MaxBoardWidth 
25 Avaiimdth=Avaiimdth-MaxBoardmdth 
} 

IFAvaiimdth>MinBoardWidth 
{ 

Select a Panel having widih=AvailWidth 
30 Avaiimdth=0 
} 

} 

} 

35 wherein: 

PositionX represents the coordinate value along first reference axis X 
shown in Figs. 4 and 5b; 

XMinSubsecb'on designates the minimum coordinate along first reference 
axis X; 

40 MaxBoardWidth designates the maximum v\ndth value that can be assigned 

to a selected panel; 
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MinBoanJWidth designates the minimum width value that can be assigned to 
a selected panel; and 

AvailWidth designates a width value associated with available surface area of 
the current subsection. 
5 Having selected the anrangement of parts providing a highest yield value, a 

basic yield value is estimated at a following step 93 represented on Fig. 1a, which 
basic yield value is associated with the arrangement of subdivided piece surface 
sections to each of which corresponds an arrangement of parts having a highest 
yield value. For so doing, all highest yield values conresponding to the anrangement 

10 of parts associated with all subsections for each piece surface section are 
cumulated. Then, whenever there is a new arrangement of subdivided piece surface 
sections to be processed, steps 65. 71, 87, 89, 91 and 93 are repeated as indicated 
by connection A-2 in Fig. la, following an affinmative result coming from condition 
checking block 102 shown in Rg. la. Whenever there is no new arrangement of 

15 subdivided piece surface sections to process, all basic yield values so obtained are 
compared to one another in order to select the arrangement of sutxlivided piece 
surface sections associated with the arrangements of parts providing a maximal 
basic yield value to be included in the optimized layout of selected parts to be cut at 
a step 85 follovwng a negative result from condition checking block 102. 

20 Turning now to Fig. lb, the optimization method according to the invention 

preferably includes a further series of steps providing improved optimization 
effidency by the consideration of a minimum yield reference value or threshold 
associated with the selected grade value i, according to whic^ the optimization 
routine explained before with reference to Rg. 1a is repeated considering at least 

25 one further optimization grade value of decreasing priority compared to the currently 
processed grade value. The minimum yield reference value may be selected 
according to the average yield expected from the grade mix of pieces of lumber and 
cun-ent optimization mode. Examples of expected average yields per standard NHLA 
grade of pieces of lumber and selected optimisation mode are given in Table 7. 

30 



Grade NHLA 


rapfirst 


Crosscut 


Flooring 


1C 


65 


63 


70 


2C 


60 


58 


67 
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3C 


55 


53 


65 



Expected yields per NHLA grade and 
optimisation mode (average %) 



TABLE 7 

5 

As shown in Rg- 1b. the further series of steps begins with a comparing step 
represented by condition checking block 104 according to which the maximal basic 
yield value is compared with a predetermined minimum yield reference value or 
threshold associated with the predetermined grade of value /, to preferably apply in 

10 the affinmative a further test represented by condition checking block 106 according 
to which the optimization process is allowed to continue as long as a predetermined 
processing time limit is not reached, in which case the optimization process is 
terminated and the obtained maximal basic yieU values are compared in a manner 
that will be later explained in detail. In other words, the basic optimization process 

15 is allowed to operate within a predetermined processing time limit until the 
arrangement of subdivided piece surface sections providing a highest maximal basic 
yield is obtained at a following step 108. Whenever the time limit condition is not yet 
satisfied, the grade value i is Incremented at step 110 prior to the repetition of step 
63, to define a further subset from the set of parts characterized by the furttier 

20 predetemnined grade value M+f associated with further predetermined group of first 
dimension values defined with respe<* to the first axis Yas shown in Fig. 3 for rip- 
first optimization mode (group of width values) or axis Xas shown in Rg. 4 for 
crosscut-first optimization mode (group of length values). Then, steps 65, 71 , 87, 
89, 91 , 93, 102 and 95 are performed again to estimate a furUier maximal yield value 

25 assodated with the fiirUier predetenmined grade value M and connesponding 
anrangement of subdivided piece surface sections, whenever the condition set forth 
at 104 is satisfied, as long as the time limit imposed at condition checking block 106 
is not reached, Othenvise, the maximal basic yield values obtained are compared 
with one another at 108 to select the arrangement of subdivided piece surface 

30 sections providing a highest maximal basic yield value to- be included in the 
optimized layout of selected parts to be cut. In other words, in a case where group 
of n predetermined grade values of decreasing cutting priority order with z=l,n , 

following a basic yield value obtained for a first predetermined grade value i=l , for 
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each predetemiined grade value i-2/i , the basic optimization sequence is repeated 
providing the condition set forth at 104 is satisfied within the time limit of condition 
106. The choice of the minimum yield reference value has a significant effect on 
grade placement sequence. A low reference value will cause an optimization 
5 process to temiinate while considering only the main selected grade, whereas a high 
reference value will cause the optimization process to consider other lower quality 
grades to increase the yield, thereby limiting placement of parts associated with the 
selected main grade. 

Whenever trie method considers more than one basic optimization mode, 

10 according an affimiatlve output generated at condition checking block 112, the 
optimization process is repeated from the initial step 62 whereby a new optimization 
mode is set as the cunrent optimization mode as indicated by connection A-l. Then, 
the following steps 63, 65, 71. 87, 89, 91, 93. 102, and 95 shown in Fig. 1a as well 
as steps 104, 106 and 108 shown in Rg. lb are repeated to obtain the 

15 conresponding arrangement of subdivided piece surface sections assodated with 
corresponding arrangement of parts providing a maximal basic y\e\6 value according 
to the selected optimization mode. Whenever there is no further optimization mode 
to consider, at a follovwng step 114, the maximum yield values of step 108 obtained 
according to all optimization modes are compared with one another, to select the 

20 arrangement of parts providing a ma)dmal basic yield value to be included in the 
optimized layout of selected parts to be cut. 

The optimization . method according to the present invention preferably 
performs a secondary part selection called "recyding" that will now be explained with 
reference to Rgs. 1c, 1d and 2, The primary part selection as explained before 

25 leaves unused surface areas within the sections as part of each piece of luml>er 
being processed. These unused surface areas may be defect adjacent area, part 
adjacent areas, or free areas that cannot be used by parts of other categories which 
cannot be considered by the primary selection due to particular defects distribution. 
The secondary selection according to the invention consist of scanning free areas 

30 remaining to the surface of a given piece section in an attempt to select parts from 
all predetermined grades, beginning with grades of higher priority. Referring now to 
Fig. 1c, a first step 116 provides the setting of a current optimization mode, 
preferably among several optimization modes to each of which is assigned a relative 
priority indicator. At step 1 1 8, a further subset of the set of parts characterized by the 
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current grade value i with i=l^n is defined. Such further subset building is required 
by the secondary selection function since the currently processed grade value is 
caused to vary. At a following step 120, there is provided data defining, for each 
arrangement of parts in each subsection included in each piece surface section 

5 providing a maximal basic yield value as explained before, one or more free areas 
according to the geometric and defect-related characteristics of the piece. The 
selected parts as part of the current optimization solution being grouped for each 
section processed, it is possible for the computer software to scan the list of parts for 
each con^sponding section to provide efficient searching. In production, the 

10 additional parts can be obtained through secondary and tertiary cutting operations 
witfi respect to first and second reference axis Y, X, for rip-first optimization mode, 
and X, Yfor crosscut-first optimization mode. 

Referring now to Rg. 15, the surface sunDunding a given selected part 
represented at 5 as been subdivided into nine areas which are grouped as follows: 

15 part pre-located areas (1, 2, 3), part located areas (4, 6) and part post-located areas 
(7, 8, 9) as part of a currently processed section 67. The general secondary part 
selection process can be implemented in the computer software using the following 
algorithm: 



20 PartStarandex^FindPartStarand€x(Order) 
PartEndIndex^FindPartEndIndex(Order) 
MinLength=FindMinLength(Order) 
MinWidih^FindMinWidth(Order) 

25 IF PartSiarandexo -1 
{ 

Sort parts for the current Solution 
NumPart=^GetNumPart(Solution) 
NtmSection-GetNtimSectionfSolution) 

30 

FROM 1=0 to i^NumSection 
{ 

Section^GetSeciion(i) 
PreceedPart-=NVLL 

35 

FROMh=0 to k=NumPart 
{ 

Part^GetPart(Solution,k) 

40 IF GetSection(Part)=^SecHon 

{ 
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IF PreceedPart=NULL 
{ 

FreeArea=PrePartArea 
Optimize FreeArea 

5 ; 

IF NOT 
{ 

FreeArea^InterPartArea 
Optimize FreeArea 

10 { 

IF Part is located near the lower edge of the section 

{ 

FreeArea-QverPartArea 
Optimize FreeArea 

15 ; 

EFNOT 

FreeArea-UnderPartArea 
Optimize FreeArea 

} 

20 PreceedPart=Part 

} 

} 

IF Part ONULL 
25 { 

FreeArea^PostPartArea 
Optimize FreeArea 

} 

IF NOT 
30 { 

FreeArea=Section 
Optimize FreeArea 

} 

} 

35 ; 



wherein: 

FindPartStartlndex(Order) designates a function allowing to retrieve the value 
of PartStartlndex according to the current cutting bill; 
40 FlndPartEndlndGx(Order) designates a function allowing to retrieve the value 

of PartEndlndex corresponding to the current cutting bill; 

FindMinLength(Order) designates a function allowing to retrieve the value of 
MinLength according to the current cutting bill; 

MInWidth represents the current minimum width value for a part; 
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FindMinWiclth(Order) designates a function allowing to retrieve the minimum 
width value for the parts included in a current cutting bill; 
NumPart represents the cunrent number of parts; 

GetNumPart(Solution) designates a function allowing to retrieve the number 
5 of part corresponding to a current processed solution; 

NumSection represents the current number of sections; 
GetNumSection(Solution) designates a function allowing to retrieve the 
number of sections involved in a currently processed solution; 
Section designates a currently processed section; 
10 PreceedPart is a parameter indicating whether there is a preceding part 

before a currently processed part and identifying such preceding part whenever it 
exists; 

GetPart(Solution, k) designates a function allowing to retrieve the data 
relating to a part k as part of a solution under processing; 
1 5 GetSection(Part) designates a function allowing to retrieve the data related to 

the section including a currently processed part; 

FreeArea designates a currentiy processed free area of raw material 
available for placing parts flierein; 

PrePartArea designates a finee area located before a part that has been 
20 incorporated in a cunrent solution through the primary selection process, 

InterPartArea designates a free area located between two parts incorporated 
into the cunrentiy processed solution through the primary selection process; 

OverPartArea designates a free area located over the part incorporated into 
the solution through the primary selection process; 
25 UnderPartArea designates a free area located under a part incorporated into 

the solution through the primary selection process; and 

PostPartArea designates a finee area located after a part that has been 
incorporated into the cumently process solution according to the primary selection 
process. 

30 In order to optimize FreeArea refenred to in the above algorithm, the 

computer software has to scan the conresponding free area while going through 
defects characterizing data related to the currently processed piece of lumber for 
checking those defects that intersect tiie considered free area in order to generate 
the larger defect-free area that can be defined within the free area, considering tiie 
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current grade limitations. Hence, if the resulting defect-free area is found to be larger 
than the minimum dimensions (width, length) of the parts included in the subset of 
parts that has been established according to the current grade and optimization 
mode, additional parts may be selected to be included within such defect-free area. 
5 The Optimize Free^rea subroutine may be programmed according to the following 
algorithm: 



Piece-GeiCurrmiFiece 
NumDefect^GetNiimDefect(Piece) 

10 

FROM i=0 to NumD^ect 
{ 

Curr€ntDefecP=GetD€fect(Piece, i) 

DefectCaiegory=GetDefeciCkitegory(Curr€ntDefect, SelectGrade, SelectFace) 
1 5 IF DefectCategory-Unacceptable OR SingleFaceAcceptable 

{ 

Extract the largest defect-free area 

} 

20 IF the defect-free area does no intersect any existing selected parts 

{ 

Place additional selected parts within the d^ect-free area 

} 

} 

25 wherein: 

Piece designates a specific piece of lumber, 

GetCurrentPiece designates a function allov^ng to retrieve identification data 
related to a currently processed piece of lumber; 

GetDefect(Piece,i) designates a function allowing to retrieve data identifying 
30 a currently processed defect /; 

Unacceptable is a parameter indicating that a given defect is considered as 
belonging to an unacceptable category of defect; 

SingleFaceAcceptable is a parameter indicating that a given defect t)elongs 
to a category according to which a defect appearing on a single face of the lumber 
35 piece is acceptable; 

The largest defect-free area is extracted in a manner that vwll be now 
explained with reference to Rg. 15. The approach used consists of identifying and 
isolating from the considered piece area, each unacceptable or single-face 
acceptable defects, for then detemnining dimensions and locations of rectangular 
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defect-free areas that maximize placement areas. It can be seen in Fig. 15 that, 
surrounding a previously selected part 5, pre-part areas includes areas 1, 2 and 3, 
over-part area includes area 4, under-part includes area 6 and post-part areas 
includes areas 7, 8 and 9. It must be pointed out that inter-part area may also be 
5 identified whenever two or more selected parts are present on a given section of a 
piece of lumber. Each of the above-mentioned group of areas will be subjected to a 
surface area calculation to establish spatial coordinates of a resulting largest defect- 
free area whose dimensions (width, length) must be larger than the minimal 
dimensions of any additional part that may be included. Once the largest defect-free 
10 area has been established, it is then possible to place one or more successive 
additional parts either horizontally or vertically according to the following algorithm: 



RepeatX^O 
RepeatY=0 

15 

AvailLength=GetLength(DrfectFreeArea) 

AvailW%dth^GetWidth(DefectFreeArea) 

NumPart=G€tNumPart(Order) 

20 FROM i^PartStartlndex to i^PartEndlndex 
{ 

Part=GetPart, (Order J) 

PartLength=GetLength(Part) 

PartWtdth^GetWidth(Part) 

25 

IFRepeatXOOamiRepeatYO 0 
{ 

30 RemainWidth=-AvailWtdth 
RemainLength='AvailLengfk 

WHILE RemainLmgth>=PartLength 
{ 

35 WHILE Remainmdth>=Partmdth 

{ 

Add a subsection 

Add a part 

} 

40 ; 

; 

; 
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wherein: 

RepeatX ls a parameter whose value corresponds to the number of times the 
, length of a current part can fit within the available length along second axis X by a 
5 currently processed defect-free area; 

RepeatY is a parameter whose value represents the number of times the 
width of a cunrent part fits within the length defined along the first reference axis Yby 
the currently processed defect-free area; 

DefectFrBeAraa is a parameter identififlng the currently processed defed-free 

10 area; 

GetLength(DefectFr&eArea) designates a function allowing to retrieve the 
length value of a currently process defect-free area; 

GetWiclth(DefectFreeArea) designates a function allowing to retrieve the 
width value of a cunrently processed defect-free area; 
15 GetNumPart(Orcler) designates a function allowing to retrieve the number of 

parts included in a cunrent cutting bill; 

RemainWidth represents the width value defined by a remaining surface 
portion of the defect-free area under processing; and 

RemainLength represents the length value defined by the remaining portion 
20 of the defe<*-free area. 

The algorithm set forth above corresponds to method step 122 shown 
in the process flow diagram of Fig. 1c, which consists of generating data defining, 
for at least some of the free areas, one or more additional parts to be included 
therein and selected from the further subset of parts, which step 122 is preferably 
25 perfonned according to a series of three sut>-steps 124, 125 and 126 executed in 
sequence as shown in Fig. 2. At first sub-step 124, data defining annangements of 
parts providing a maximal basic yield value is translated in a first direction along the 
second reference axis within conresponding subsections to extend the free areas, 
and the additional parts from the further subsets of parts are selected accordingly. 
30 The effect of such translation in the context of a rip-first optimization mode is shown 
in Fig. 16 as part of a computer display screen, wherein each of selected parts 128 
and 129 as a result of the primary selection process has been translated in the 
direction of axis in order to extend free areas 130, 131 respecHvely adjacent 
selected parts 128, 129, prior to proceed with the selection of the additional parts 
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from the further subset of parts as described before. In the example shown in Fig. 
16, none of the listed remaining parts has been found to satisfy the dimensional 
requirements characterizing defect-free areas 130 and 131, and therefore the 
secondary selection process does not contribute to the total basic surface yield 
5 obtained through the primary selection. At following sub-step 125 shown in Fig. 2, 
the resulting data is further translated in a second direction opposed to the first 
direction along the second reference X axis within corresponding subsections to 
further extend the free areas, prior to selecting new additional parts from the further 
subset of parts as explained before. As shown in Fig. 17, the translation of selected 

1 0 parts 1 28. 1 29 in the left direction toward the origin of axis X gives the opportunity to 
use further defect-free areas 132, 133 respectively adjacent to selected parts 128, 
129. Sometimes, selected parts translation to the left may not necessarily provide 
surface yield increase when none of the remaining parts satisfies the dimensional 
requirements characterized by any of further defect-free areas 132 or 133, However, 

15 any insertion of additional parts resulting from the first translation step 124 may have 
an effect on following translation steps, and the intermediate translation step 125 is 
prerequisite to the following translation step 126 which will be now explained with 
reference to Fig. 18. At step 126 shown in Fig. 2, further translating the data 
defining one ore more of said parts and resulting from said step b) in said first 

20 direction within conresponding said subsections to further extend said free areas, 
and selecting new said additional parts firom said further subset of parts accordingly. 
The data defining part 129 resulting from prior translation step 125, is further 
translated in the first direction, i.e. to the right along axis AT. within conresponding 
subsection, allowing to selecting new additional parts 135, 136 from the further 

25 subset of parts to be included within the defect-free area extended as a result of the 
further translation applied at step 126, as new way to further extend available inter- 
subsection surface areas. It Is to be understood that the further translation of step 
126 in the first direction may be applied to more than one part, dependent on the 
number of subsections involved. Rg. 19 represents a cutting solution display screen 

30 generated by the computer software in the context of another optimization process 
example, showing the resulting selected parts layout for both faces 117, 117' along 
with specific infomiation relate to each selected part in the form of a table 119. The 
computer software may also be program to display other complementary data 
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related to specific selected part, such as production identification number, position 
data (Xmin . Xmax , Ymin . Ymax ) and surface data. 

Turning now to Figs. 7 and 8, the two-face approach for the segmentation of 
each piece surface section into subsections will be now explained in grater detail. 
5 In order to perfonn accurate generation of anangements of parts to be included in 
each subsecb'on generated by the process, while incorporating into the layout 
highest grade parts, an analysis of the locations of deifects is performed using the 
algorithm set out above, which considers both faces 73, 73' of a piece of lumber 
under processing. For so doing, a function according to which defects associated 

10 with a particular category known as single-face category can be considered by the 
computer software as desired by the operator.. As can be seen on Figs. 7 and 8, 
each piece section of raw material shown has first and second opposed faces 73, 
73' which are preferably chosen as the main faces of the piece of lumber, wherein 
the data representing geometric and defect-related characteristics of the piece 

15 generated in step 71 of the flow diagram of Fig. la, includes data characterizing 
each effaces 73. 73' wherein defect-related characteristics are defined according to 
predetemnined defect types, each of which is associated with one of a plurality of 
defect acceptability categories including acceptable defect, single-face acceptable 
defect and an unacceptable defect categories as listed in the table displayed on the 

20 computer screen shown in Fig. 6, as explained before. In this context, the data 
generating step 71 shown in Rg. la is preferably performed so that substantially 
none of the subsections designated at 81 , 83, 85 for first face 73 and by 81', 83', 85' 
for second face 73' of the piece section represented in Fig. 7, includes two (2) 
defects associated with the single-face category respectively present on first and 

25 second faces 73, 73', while sharing a same portion or subsection 83 of the piece 
section of raw material. It can be seen from Fig. 7 that none of the single-face 
acceptable defects 96 appearing onto first face 73 has corresponding defects 
associated with the single-face category present onto the second face 73' as part of 
the same subsection designated at 83, 83'. It can also be seen that the single-face 

30 acceptable defect 96' appearing on second face 73' does not correspond to any 
such single-face acceptable defect appearing onto first face 73 at the same 
subsection designated at 85, 85'. 

Tuming now to Fig. 8, it can be seen that the above-mentioned condition is 
not satisfied in respect of single-face acceptable defects designated at 97, 97 
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appearing respectively onto first face 73 and second face 73' of the piece section 
shown, since the two single-face acceptable defects 97, 97* share a same section 
area designated at 88, 88'. The computer software causes the shared portion to be 
removed by placing a secondary cutting axis 101 at 

Xmax location with respect 

5 to AT axis shown in Fig. 8. It can be appreciated that the remaining defect portions 
99, 99' respectively assodalsd with subsections 83 and 85* now comply with the 
above-mentioned condition. It can t>e seen from both Figs. 7 and 8 that the 
presence of an unacceptable defect 100 on first face 73 causes the computer 
software to generate secondary cutting lines 75, 77 even if no defect appears onto a 

10 corresponding portion 102 of second face 73. All unacceptable defects 
systematically cause the computer software to generate a secondary cxjtting line to 
define the end portion of a current portion, and a further secondary cutting line to 
define the beginning of a following subsection. However, vyrfien a first encountered 
single-face acceptable defect is identified as well as the specific face onto such first 

15 defect is found, the computer software does not generate a secondary cutting line. 
Any following single-face acceptable defects appearing on the same face as that of 
first find single-face acceptable defect are ignored until another single-face 
acceptable defect is found onto the opposed face so as to generate a secondary 
cutting line to complete the current subsection. In other words, in cases where two or 

20 more defects are aligned one with another while being not present onto the same 
face, as shown in the example of Fig. 8, the Xmax 'imft assodated with the defect 
present on a given face which precedes an aligned defect present onto the other 
face is used to detemiine the next border of the current subsection- 

Tuming back to Rg. 1c, following data generating step 122 is a 

25 comparing step represented by condition checking block 135 according to which the 
current grade value /, is compared to the predetemriined numt>er n of grades, and 
whenever the such limit is not reached, the current grade value i is incremented at 
step 137 prior to the repetition of step 118 through connection C-2, to define a 
further subset from the set of parts characterized by the further predetermined grade 

30 value /=/+^. Data generating steps 120 and 122 are then performed in a same 
manner as explained before, and the condition set fortti in block 135 is applied again 
with the current value of /. Whenever /=/?, additional yield values associated with the 
inclusion of the additional parts resulting from the recycling function applied to the 
cunrent optimization mode for the n grade values are estimated at step 139. Tuming 
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now to Fig. 1d, whenever the recycling function considers more than one basic 
optimization mode, according an affirmative output generated at condition checking 
block 141, the sequence of recycling steps is repeated from the initial step 116 
whereby a new optimization mode is set as the cunBnt optimization mode as 
5 indicated by connection C-1. Whenever there is no further optimization mode to 
consider as tested at 141, resulting additional yield values associated with all tested 
optimization modes are compared at step 143, to select the additional parts 
providing a maximal additional yield value. 

Referring now to Fig. ie, an optional function allowing the 

10 incorporation, into the optimized solution providing a maximal yield value, of furUner 
additional elements of various dimensions, which function is called a reclaiming 
function, will now be presented- Primary placement of parts as well as seopndary 
placement of recycled parts may still leave free areas made of unused raw material, 
within some sections of a processed piece of lumber. The reclaiming function 

15 programmed in the computer software consists of inserting one or more additional 
elements of various dimensions, within predetermined limits, whenever the criteria 
associated witii the reclaimed grade, i.e. whether particular types of defects present 
onto the section under processing are acceptable, are met. Since predetermined 
dimensional characteristics are not specified, the set of parts defined by the cunrent 

20 cutting bill need not to be considered. In otiier words, the reclaiming function 
involves an alWehgth/all-width approach for additional element placement. More 
specifically, as shown in Fig. 1e, the reclaiming function involves a step 145 
consisting of generating data defining, for each arrangement of parts in each 
subsection included in each piece surface section and providing a maximal basic 

25 yield value, one or more remaining free areas according to the geometric and defect- 
related characteristics of the piece. Then, ttiere is provided a step 147 consisting of 
generating data defining, for at least some of the remaining free areas, one or more 
reclaimed elements to be included therein as part of the optimized layout of selected 
parts to be cut Finally, additional yield values associated with the inclusion of tiie 

30 reclaimed elements can be estimated at step 149. 

A mentioned before, the computer software has access to a replacement 
parts data file selected for the cunent simulation processing, which is used to 
replace a part having been selected in a layout, and therefore vwthdrawn from the 
active set of parts, by a new required part listed in the current cutting bill, taken from 
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an inactive portion thereof. Such replacement function may be perfonmed according 
to the replacement priority indicator value given to each part type, on the basis of 
anyone of the following rules: (1) a replacing part must be of a same grade, same 
production and of substantially a same dimension (length or width) vwthin a 
5 predetennined dimension range, than those of the replaced part; (2) a replacing part 
must be of a same grade, same production, and of the nearest dimension (length or 
v^dth), than those of the replaced part; (3) a replacing part must be of a same grade 
and of substantially a same dimension (length or width) within a predetenmined 
dimension range, than those of the replaced part; and (4) a replacing part must be of 

10 a same grade and of the nearest dimension (length or width), than those of the 
replaced part. It must be understood other rules allowing that a replacing part be of 
higher or lower grade, or allowing a replacing part be of a different type, may be 
established. Preferably, the replacement function must ensure to maintain a 
minimum, predetermined number of parts of each type in the active set One 

15 approach to meet such object consists of applying weighing factors to the estimated 
area of each selected part, whose values are chosen within a "0 - 1" range 
depending upon the priority given to each type of parts. In the particular case of 
panel cutting optimization, a high weighing factor value for paneling may Increase 
the resulting yield while causing the production of fixed-sized components to 

20 decrease. In such case, a lower value for the minimum yield reference would allow 
faster placement of components. Weighing factors associated with indrvidual parts 
can also be considered in the comparison of optimization trial solutions, wherein the 
maximal basic yield values thereof are compared one with another to select the 
arrangement of subdivided piece surface sections providing a highest maximal basic 

25 yield value, as explsuned before. The weighing factors are used along with 
comparison rules according to which the primary and secondaiy (recycled) yield 
values, or a combination thereof, is considered in the selection of the final solution, 
either on the basis of surface or economic value yields. 

Turning now to Fig. 20, a simulation results display screen generated by 

30 computer software is shown, giving a list of parameters related the cunent 
simulation process, including selected grade for primary optimization and resulting 
yield data for each grade considered. Refem'ng to Fig. 21, a production results 
display screen generated by the computer software is shown, giving simulation 
related data such as name of simulation, optimization mode, number of pieces 



wo 2004/051523 



PCT/CA2003/001872 



47 



involved, number of parts selected, input surface of pieces, part surface generated 
and monetary value ($), along with a table presenting production data related to 
each part selected, namely: part code, product type, quantity produced, surface 
generated, surface percentage of input surface, and piece input surface. 
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We claim : 

1 . A method of optimizing a layout of selected parts to be cut from a piece of raw 
material with respect to first and second orthogonal reference axis, using data 
representing geometric and defect-related characteristics of said piece and data 
representing geometric and grade characteristics of a set of parts to be cut, said 
method comprising the steps of: 

i) defining a subset of said set of parts characterized by a predetermined 
grade value and assodated with a predetermined group of first dimension values 
defined with respect to said first axis; 

ii) generating data defining an arrangement of sulxlivided piece surface 
sections to be obtained through a primary cutting operation with respect to said 
second reference axis and according to one or more of said first dimension values; 

iii) generating data defining one or more subsections included in each said 
piece surface section according to said geometric and defect-related characteristics of 
said piece; 

iv) generating data defining for each said subsection a plurality of 
arrangements of parts to be included therein and selected firom said subset of parts, 
to be obtained through a secondary cutting operation with respect to said first 
reference axis; 

v) estimating yield values associated with said arrangements of parts; 

vi) comparing said yield values to select one of said annangements of parts 
having a highest yield value; 

vii) estimating a basic yield value for said arrangement of subdivided piece 
surface sections; 

viii) repeating said steps ii) to vii) for one or more new amangements of 
subdivided piece surface sections to estimate corresponding basic yield values; and 

ix) comparing all said basic yield values one with another to select said 
arrangement of subdivided piece surface sections assodated with said arrangements 
of parts providing a maximal basic yield value to be included in said optimized layout 
of selected parts to be cut 

2. The method of daim 1, wherein said piece is an elongate piece having its 
longitudinal axis parallel to said second reference axis, said piece having its 
transverse axis parallel to said first reference axis. 
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3. The method of daim 1, wherein said piece is an elongate piece having its 
longitudinal axis parallel to said first reference axis, said piece having its transverse 
axis parallel to said second reference axis. 

4. The metiiod of claim 1 , further comprising the steps of: 

x) comparing said maximal basic yield value with a predetermined 
minimum yield reference value associated with said predetermined grade, and, 
whenever said maximal basic yield value is lower than said minimum yield reference 
value, performing the steps of: 

xi) defining a further subset of said set of parts characterized by a further 
predetenDined grade value and associated with a further predetermined group of first 
dimension values defined with respect to said first axis; 

xii) repeating said steps ii) to ix) to estimate a further maximal yield value 
associated with said further predetermined grade value for corresponding said 
arrangement of subdivided piece surface sections ; 

xiii) comparing said maximal basic yield values one with another to select 
said anrangement of subdivided piece surface sections providing a highest maximal 
basic yield value to be included in said optimized layout of selected parts to be cut. 

5. The method of claim 4, wherein said step ii) to xiii) are performed within a 
predetermined processing time limit until said anrangement of subdivided piece 
surface sections providing a highest nnaximal basic yield is obtained. 

6. The method of daim 1 , wherein said predetenmined grade value is a first one 
of a group of n predetemriined grade values of decreasing cutting priority order i=l^ , 
for each said predetermined grade value i=2^ , the method further comprising the 
steps of: 

x) comparing said maximal basic yield value with a predetermined 
minimum yield reference value associated with said first predetermined grade, and, 
whenever said maximal basic yield value is lower than said minimum yield reference 
value, performing the steps of: 

xi) defining a further subset of said set of parts characterized by said 
predetermined grade value i and associated with a corresponding further 
predetermined group of first dimension values defined with respect to said first axis; 
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)di) repeating said steps ii) to ix) to estimate a further maximal yield value 
associated with said predetermined grade value i for conresponding said 
arrangement of subdivided piece surface sections; 

xiii) comparing said maximal basic yield values one wiUi another to select 
said anrangement of subdivided piece surface sections providing a highest maximal 
t)asic yield value to be included in said optimized layout of selected parts to be cut. 

7. The method of claim 6, wherein said step ii) to xiii) are performed within a 
predetemnined processing time limit until said arrangement of subdivided piece 
surface sections providing a highest maximal basic yield is obtained. 

8. The method of claim 1, wherein said piece is an elongate piece defining 
longitudinal and transverse axis, said first and second orUiogonal reference axis 
defining one of a first and second optimization mode, said longitudinal and transverse 
axis being respectively parallel to said second and first reference axis according to 
said first optimization mode, said longitudinal and transverse axis being respectively 
parallel to said second and first reference axis according to said second optimization 
mode, said steps i) to ix) being first perfonned according to said first optimization 
mode, the metiiod furtiier comprises the steps of: 

x) selecting said second optimization mode; 

xi) repeating said steps i) to ix) to obtain corresponding said arrangement 
of subdivided piece surface sections associated with corresponding said 
arrangements of parts providing a maximal basic yield value according to said second 
optimization mode; and 

xii) comparing the maximum basic yield values resulting from said steps 
ix) and xi) obtained according to said first and second optimization modes 
respectively, to select said arrangement of parts providing a maximal basic yield value 
to be included in said optimized layout of selected parts to be cut. 

9. The method of daim 1, wherein said predetennriined grade value is a first one 
of a group of n predetenmined grade values of decreasing cutting priority order i^l^t , 
for each said predetennined grade value i=M . the method further comprising tiie 
steps of: 

x) defining a furtiier subset of said set of parts characterized by said 
predetenr^ined grade value / ; 
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xi) generating data defining, for each said arrangement of parts In each 
said subsection included in each said piece surface section and providing a ma)dmal 
basic yield value, one or more free areas according to said geometric and defect- 
related characteristics of said piece; 

xii) generating data defining, for at least some of said free areas, one or 
more additional parts to be included therein and selected from said further subset of 
parts, as part of said optimized layout of selected parts to be cut, 

1 0. The method of daim 9, further comprising the step of: 

xiii) estimating additional yield values associated with the inclusion of said 
additional parts. 

11. The method of claim 9, wherein said additional parts are to be obtained 
ttirough secondary and tertiary cutting operations witii respect to said first and second 
reference axis, respectively. 

12. The mettiod of daim 9, wherein said piece is an elongate piece having its 
longitudinal axis being parallel to said second reference axis, said step xii) induding 
the steps of : 

a) translating said data defining said anrangements of parts providing a 
maximal basic yield value in a first direction along said second reference axis within 
corresponding said subsections to extend said free areas, and selecting said 
additional parts from said further subset of parts accordingly; 

b) translating tiie data resulting from said step a) in a second direction 
opposed to said first direction along said second reference axis within corresponding 
said subsections to further extend said free areas, and selecting new said additional 
parts from said furUier subset of parts accordingly; and 

c) furttier translating the data defining one ore more of said parts and 
resulting from said step b) in said first direction within corresponding said subsections 
to further extend said free areas, and selecting new said additional parts from said 
furtiier subset of parts accordingly. 

13. The method of daim 12, wherein said additional parts are to be obtained 
through secondary and tertiary cutting operations witti respect to said first and second 
reference axis, respectively. 



wo 2004/051523 



52 



PCT/CA2003/001872 



14. The method of claim 9. further comprising the steps of: 

xiii) generating data defining, for each said arrangement of parts in each 
said subsection included in each said piece surface section and providing a maximal 
basic yield value, one or more remaining free areas according to said geometric and 
defect-related characteristics of said piece; and 

xiv) generating data defining, for at least some of said remaining free 
areas, one or more reclaimed elements to be included therein as part of said 
optimized layout of selected parts to be cut. 

15. The method of daim 14, further comprising the step of: 

XV) estimating additional yield values associated with the inclusion of said 
reclaimed elements. 

16. The method of daim 9. wherein said predetermined grade value is a first one 
of a group of n predetemnined grade values of decreasing cutting priority order , 
wherein said piece is an elongate piece defining longitudinal and transverse axis, said 
first and second orthogonal reference axis defining one of a first and second 
optimization mode, said longitudinal and transverse axis being respectively parallel to 
said second and first reference axis according to said first optimization mode, said 
longitudinal and transverse axis being respectively parallel to said second and first 
reference axis according to said second optimization mode, said steps i) to ix) being 
first performed according to said first optimization mode, the method further 
comprising, for each said predetermined grade value , the steps of: 

x) defining a further subset of said set of parts characterized by said 
predetermined grade value i ; 

)d) generating data defining, for each said arrangement of parts in each 
said subsection induded in each said piece surface section and providing a maximal 
basic yield value, one or more free areas according to said geometric and defect- 
related characteristics of said piece; and 

xii) generating data defining, for at least some of said free areas, one or 
more additional parts to be included therein and selected from said furtiier subset of 
parts; said method furtiier comprising the steps of: 

xiii) estimating additional yield values assodated with said additional parts; 

xiv) repeating said steps x) to xii) according to said second optimization 
mode and estimating additional yield values associated witii corresponding said 
additional parts; 
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XV) comparing the additional yield values resulting from said steps xili) and 
xiv) respectively, to select said additional parts providing a maximal additional yield 
value, to be included in said optimized layout of selected parts to be cut. 

17. The method of claim 16, further comprising the steps of: 

xvi) generating data defining, for each said anrangement of parts in each 
said subsection included in each said piece surface section and providing a maximal 
basic yield value, one or more remaining free areas according to said geometric and 
defect-related characteristics of said piece; and 

xvii) generating data defining, for at least some of said remaining free 
areas, one or more reclaimed elements to be included therein as part of said 
optimized layout of selected parts to be cut. 

1 8. The method of claim 1, wherein said piece of raw material has at least first and 
second opposed faces, said data representing geometric and defect-related 
characteristics of said piece including data characterizing each said first and second 
opposed faces, said defect-related characteristics being defined according to 
predetermined defect types each of v\^ich being associated with one of a plurality of 
defect acceptability categories including acceptable defect, single-face acceptable 
defect and unacceptable defect categories, said data generating step iii) being 
performed so that substantially none of said subsections includes two defects 
assodated with the single-face category and respectively present on said first and 
second faces. 

19. A software product data recording medium in v^ich program code is stored, 
which program code will cause a computer to perform a method of optimizing a layout 
of selected parts to be cut from a piece of raw material with respect to first and 
second orthogonal reference axis, using data representing geometric and defect- 
related characteristics of said piece and data representing geometric and grade 
characteristics of a set of parts to be cut, said method comprising the steps of: 

i) defining a subset of said set of parts characterized by a predetermined 
grade value and associated with a predetermined group of first dimension values 
defined with respect to said first axis; 

ii) generating data defining an arrangement of subdivided piece surface 
sections to be obtained through a primary cutting operation with respect to said 
second reference axis and according to one or more of said first dimension values; 
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iii) generating data defining one or more subsections included in each said 
piece surface section according to said geometric and defect-related characteristics of 
said piece; 

iv) generating data defining for each said subsection a plurality of 
arrangements of parts to be included therein and selected from said subset of parts, 
to be obtained through a secondary cutting operation with respect to said first 
reference axis; 

v) estimating yield values associated with said arrangements of parts; 

vi) comparing said yield values to select one of said arrangements of parts 
having a highest yield value; 

vii) estimating a basic yield value for said arrangement of subdivided piece 
surface sections; 

viii) repeating said steps ii) to vii) for one or more new arrangements of 
subdivided piece surface sections to estimate conresponding basic yield values; and 

ix) comparing all said basic yield values one with another to select said 
arrangement of subdivided piece surface sections associated with said anrangements 
of parts providing a maximal basic yield value to be included in said optimized layout 
of selected parts to be cut. 
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